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Abstract 

Time variability of the solar neutrino flux especially in the low and intermediate 
energy sector remains an open question and, if it exists, it is likely to be originated 
from the magnetic moment transition from active to light sterile neutrinos at times 
of intense solar activity and magnetic field. We examine the prospects for the SNO+ 
experiment to address this important issue and to distinguish between the two classes 
of solar models which are currently identified as corresponding to a high (SSM I) 
and a low (SSM II) heavy element abundance. We also evaluate the predictions 
from these two models for the Chlorine experiment event rate in the standard LMA 
and LMA+Resonant Spin Flavour Precession (RSFP) scenarios. It is found that 
after three years of SNO+ data taking, the pep flux measurement will be able to 
discriminate between the standard LMA and LMA+RSFP scenarios, independently 
of which is the correct solar model. If the LMA rate is measured, RSFP with 
Bq ~ 280kG for the resonant Am^i can be excluded at more than 4c7. A low rate 
would signal new physics, excluding all the 90% allowed range of the standard LMA 
solution at 3cr, and a time variability would be a strong signature of the RSFP 
model. The CNO fluxes are the ones for which the two SSM predictions exhibit 
the largest differences, so their measurement at SNO+ will be important to favour 
one or the other. The distinction will be clearer after LMA or RSFP are confirmed 
with pep, but still, a CNO measurement at the level of SSM I/LMA will disfavour 
SSM II at about 3 a. We conclude that consistency between future pep and CNO flux 
measurements at SNO+ and Chlorine would either favour an LMA+RSFP scenario 
or favour SSM II over SSM I. 
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1 Introduction 



Neutrino oscillations in matter [Ij with the resonant amplification of a small vacuum 
mixing angle [2j, although a much attractive mechanism, has not proven to be the origin 
of the solar neutrino deficit. While also an oscillation and resonant effect, the large mixing 
angle solution (LMA) [31 HI [5] has instead become generally accepted as the dominant one 
[6l[71[8]. In the LMA mechanism the conversion from active electron neutrinos produced in 
the solar core to weakly interacting ones of another flavour takes place through a strongly 
adiabatic resonance occuring still at the solar core. The order parameter is a large vacuum 
mixing of the order of 30° — 33°. 

On the other hand, the time variability of the active neutrino event rate had been 
hinted long ago by the Chlorine collaboration p] who suggested a possible anticorrelation 
of active neutrino flux with sunspot activity. It was then interpreted by Voloshin, Vysotskii 
and Okun pil] as a neutrino magnetic moment effect, such that an intense sunspot activity 
would induce a conversion of a large fraction of neutrinos into undetectable ones (either 
steriles or of a different flavour) through the interaction of the magnetic moment with the 
solar magnetic field. Hence a more intense solar activity would correspond to a smaller 
flux of detectable neutrinos and viceversa. An interesting evolution of this proposal was 
the suggestion in 1987 by Lim and Marciano and by Akhmedov [12\ that the neutrino 
spin flavour precession could take place anywhere inside the sun via a resonant process. 
This enhances the mechanism and allows for a smaller neutrino magnetic moment in 
order to produce a visible effect. The resonant spin flavour precession (RSFP) bears a 
resemblance to matter oscillations and is a result of the balance between matter density 
and the product ^^B (neutrino magnetic moment times the solar field). 

The interpretation of solar data is at present still partly ambiguous, with several sce- 
narios involving RSFP [TSj [H] and non-standard neutrino interactions [15] being viable. 
Our knowledge of the solar neutrinos relies essentially on the data from the high energy 
sector (mainly the flux), whereas the overwhelming low and intermediate energy one 
remains vastly unknown, except for the integrated measurements provided by the radio- 
chemical experiments. As pointed out earlier [TB], the observed decrease of the Gallium 
event rate |T7j opens the question of whether there is time variability affecting only the 
low energy sector and has motivated the investigation of RSFP to light sterile neutrinos in 
combination with LMA 0. LMA+RSFP thus requires a sizable neutrino magnetic moment 
and a strong field at times of intense solar activity and a weak field otherwise, which causes 
the modulations in the neutrino event rate. Low energy solar neutrino experiments like 
Borexino [TSl US], KamLAND [Bl [7] and SNO+ will no doubt help in clarifying the situ- 
ation within the next few years. The first Borexino data, very recently released [19], are 
compatible with the RSFP model predictions [13l [16] earlier derived since those data were 
taken during the present year (2007) when the magnetic solar activity is at a minimum. 

Besides our limited knowledge of the low energy neutrino sector, one of the key inputs 

^The magnetic moment we refer here is of course a Majorana transition one since it connects active 
neutrinos to sterile ones, hence of a different flavour. 
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of solar models, namely the amount of heavy element abundance relative to hydrogen, 
ZjX^ is still unclear [20]. This affects the neutrino fluxes, in particular most strongly the 
CNO ones and to a lesser extent the ^ Be one. Two classes of standard solar models (SSMs) 
may at present be distinguished [211 [22], one with a 'high' value of ZjX [23j and another 
with a 'low' one [21], hereby denoted by SSM I and SSM II, respectively. In the first the 
metallicity is consistent with sound speed, convective zone depth and density profiles, in 
excellent agreement with helioseismology. Such is not the case in the second one which 
is however based on an improved modeling of the solar atmosphere. As uncertainties are 
considerable, we will in the present paper take both models into account H. 

The purpose of this paper is to analyse the possibility of the proposed SNO+ experi- 
ment to ascertain on whether RSFP occurs at low and/or intermediate energies. We will 
also investigate the possible signatures of the two classes of SSMs, in particular whether it 
is possible to distinguish between them with data from the forthcoming SNO+ and from 
the Chlorine ^25j experiment. For neutrino-electron scattering experiments, sensitivity to 
neutrino physics depends on the accuracy of solar models, since the calculation of the 
electron neutrino survival probability relies on the comparison of the measured flux with 
the total predicted flux. So for the purpose of distinguishing models with different survival 
probabilities, the best sensitivity lies in the observation of the solar flux component with 
the smallest error. Above the threshold of liquid scintillator electron scattering experi- 
ments, the component with the most accurate prediction is the monoenergetic (1.442 MeV) 
source of pep neutrinos. The pep flux also has the smallest spread in the predictions from 
the main solar models, which allows its measurement to be sensitive to Neutrino physics 
without needing to choose a particular solar model. 

The paper is divided as follows: in Section [2] we briefly review the RSFP model and 
its motivations, and compare its predictions with standard LMA for the existing Chlorine 
experiment data, in the context of the two SSMs. In Section [3] we briefly describe the 
SNO-I- experiment. In Sec. 14.11 we describe the method and present our results. In 
Section 14.21 we comment on the sensitivity of the pep measurement to an RSFP effect 
and we investigate the possibility for SNO+ to distinguish between SSMs of the two types 
with the CNO measurement in Sec. 14.31 Finally in Section [5] we report our conclusions. 

2 Resonant spin flavour precession and solar models 

In this Section we describe the RSFP effect from active to light sterile neutrinos and 
evaluate the survival probability for LMA and LMA+RSFP (non-vanishing [i^ with low 
and high solar magnetic field). We start with a brief review of the model presented in ref. 
[TBj whose predictions for Borexino and LENS were analyzed in ^13j and for KamLAND in 
[26] . We also evaluate the event rates in the Chlorine experiment for standard neutrinos 
(i.e. massless and with no magnetic moment), for the LMA and LMA+RSFP scenario 
and for standard LMA (vanishing We use the neutrino fluxes from the type I [23] 
and II [211 solar models. 

^It should be mentioned however that SSM I is referred to as the 'preferred' solar model [2T] . 
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2.1 Resonant Spin Flavour Precession 



The possible anticorrelation with sunspot activity of the electron neutrino flux in the Ga 
experiments is most naturally explained in terms of a resonant conversion to neutrinos 
of other types that are unseen by the weak charged current. As shall be discussed, this 
requires a mass square difference between the intervening neutrino flavours in the resonance 
ArTT-Q^ = O(10~^e\^^). Such a value implies that conversion to weakly interacting neutrinos 
is excluded, leaving us the possibility of conversion to sterile neutrinos. The simplest 
departure from conventional LMA able to generate such a conversion is provided by a 
model which, in addition to the two flavours involved in LMA, introduces a sterile neutrino 
with a vanishing vacuum mixing. The active states z/g, i/^ communicate to the sterile 
one via a single magnetic moment. Owing to the large order of magnitude difference 
between the parameters Am|]^ and AmQ^^ (associated with the LMA and SFP resonances 
respectively) the two resonances are located far apart, so that they do not interfere. A 
straightforward but long calculation leads to the following form of the Hamiltonian [16] 



in the mass matter basis (i>o i^i ^2)- In this equation l^, Vx are the matter induced 
potentials for z/g and z/^, B is the solar magnetic field and 9 is the vacuum mixing angle. 

The important transition with order parameter iiyB and whose time dependent effi- 
ciency may determine the possible modulation of neutrino flux is therefore between mass 
matter eigenstates z/q, ui- It is expected to resonate in the region where the magnetic fleld 
is the strongest in the period of high solar activity. We consider the Landau Zener ap- 
proximation in dealing with the two resonances. Since the LMA one is strongly adiabatic, 
we need only to consider the jump probability between z/q, z>i in the vicinity of the SFP 
resonance 

In the LMA-(-RSFP scenario active neutrinos are partially converted to light sterile 
ones at times of strong solar magnetic fleld, thus leading to the lower Gallium event 
rate in the period 1998-2003 [17J, while LMA acts alone otherwise and the higher rate 
is obtained. As in previous publications, throughout our RSFP calculations we use a 
value for the neutrino magnetic moment ^i, = 10~^^^b- For a deflnite neutrino energy 
the critical density is flxed by the order of magnitude of the corresponding mass squared 
difference and determines the resonance range. Furthermore it has been noticed [21] that 
the solar rotation frequency matches the observed neutrino modulation in the equatorial 
section of the convection zone near the tachocline, at ~ (0.7— 0.8)i?5. This indicates a fleld 
proflle whose time dependent peak occurs around this depth. In order that the low and 
intermediate energy neutrino SFP resonance is located in this range, so as to provide time 
modulation in this sector, one needs AttIq^ = O(10~^)el^^. Hence, owing to the absence of 
interference between the two resonances (LMA and RSFP) the high energy solar neutrino 

•^For calculational details see ref. [T6] . 
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Rj (SNU) 




SSM I 


SSM II 


SSM I 


SSM II 


pep 


1.42 X 10^(1 ±0.015) 


1.45 X 10^(1 ±0.01) 


0.222 ±0.003 


0.226 ±0.002 


''Be 


4.84 X 10^(1 ±0.105) 


4.34 X 10^(1 ±0.093) 


1.16 ±0.122 


1.043 ±0.097 


85 


5.69 X 10^(1 ±0.16) 


4.51 X 10^(1 ±0.12) 


6.740 ±1.078 


5.342 ±0.641 




3.05 X 108(l±0:i) 
2.31 X 10*(1±°-|) 


2.00 X 108(l±[j:i45) 
1.44 X 108(l±0;i65) 


0.052 ±H1! 
0.154 ±lf,l 


0.034 ±Hg! 
0.096 ±',tl 


^Rci 


8.33 ± 1.21 


6.74 ±0.73 



Table 1: Predicted total fluxes (0j) and expected event rates (Rj) in the Chlorine experi- 
ment for standard neutrinos. 



experiments (SuperKamiokande [28] and SNO [8]) are not expected to exhibit any time 
modulation in their event rate. On the contrary, of all the experiments so far, Gallium 
will be the most sensitive of all to variability, while Chlorine may lie in the borderline 
with some moderate variability which failed to be clearly detected. 

Based on the above criteria we chose solar fleld proflles as in flg.l of ref. [13]. In 
our previous publications [TH [131 130] we performed two separate flts for the high and low 
Gallium event rates together with all other solar neutrino data. We also investigated the 
dependence of the results on the choice of solar fleld proflle and on the value of Amg^. We 
will in the present paper use our best choices from ref. ^30j , namely 

B= ^ 0<X<Xe (2) 

cosn[b[x — Xc)\ 

B = TFTTT ^ Xc < X < 1 (3) 

cosh[lo{x — Xc)\ 



for the fleld proflle and 



Aml-^ = -1.7 X IQ-^eV^. 



Here Bq is the peak fleld value, which for active sun we take to be 280A;G, x is the fraction 
of the solar radius and Xc = 0.71. 
Furthermore we use [291 [30] 



e = 0.508, Am' = 8.2 x lO^^eV"^ (4) 

which lie within la of the central values for KamLAND[^. The fleld proflle chosen, 
together with the value of Am^^, will be responsible for a modulation shared by pp and 
^Be neutrinos and a dip at the low and intermediate energy in the survival probability. 
This is shown in flg[Il 



2.2 Chlorine rate, type I and II solar models 



4 




Figure 1: Electron neutrino survival probability for LMA (zero field) and LMA+RSFP 
for the considered point. For the RSFP case we use a peak field value Bq = 280kG and a 
profile as in [SUED]. 



Here we analyse the Chlorine event rates in the two scenarios, the one with and the 
one without neutrino magnetic moment, using the neutrino fluxes from SSM I [23] and 
SSM II |24j. It will be seen, from the comparison of these predictions with the Chlorine 
data [25] 

Rci = 2.56 ± 0.16 ± 0.15 SNU, (5) 

that RSFP is compatible with both solar models while standard LMA favours SSM II. The 
predicted fluxes (pj, the partial event rates Rj for each flux and the total event rate for the 
Chlorine experiment Rci are given in table ([T]) for both models with standard neutrinos. 
The errors in the partial event rates 6Rj listed in this table were calculated from 

Rj = J a^, (j)^ (l ± dE, = Rj ± SRj (6) 

for flux j that is, only the flux errors were considered. 

The results for the event rates are given in table ([2]) for LMA with parameters as in 
eq. (jll) and LMA+RSFP. The errors in the total Chlorine rate are obtained using the 
correlation among the flux errors as in tables 16 and 17 of ref. [21]. Comparing the event 
rates for LMA (zero field) and LMA + RSFP (high field), it is seen that some modulation 
is expected in both SSMs, due to the time dependence of the solar magnetic field. 
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Rlma (SNU) 


i?i?5FP (SNU) 




SSM I 


SSM II 


SSM I 


SSM II 


pep 
'Be 
85 

15^ 


0.133 ±0.002 
0.710 ±0.075 
2.003 ±0.320 

n noo +0.010 

U.UOZ ICq QQg 

Q nqi -uO.030 
u.uy± 3ZQQ26 


0.136 ±0.001 
0.637 ±0.059 
1.587 ±0.191 
0.021 ±0.003 

n nc;7 _(_0.009 
U.UO^ ='=0.008 


0.088 ±0.001 
0.447 ±0.047 
1.888 ±0.302 

0.019 ±0:»i 

0.059 ±rofr 


0.090 ±0.001 
0.401 ±0.037 
1.497 ±0.180 
0.012 ±0.002 

n nQ7 _(_0.006 


^Rci 


2.97±0.40 


2.44±0.25 


2.50±0.35 


2.04±0.22 



Table 2: Expected event rates in the Chlorine experiment for the two SSMs. For the 
parameter choices used in the LMA and LMA±RSFP cases, see the main text [eq. (jl])] 
and fig. [H 



Averaging the results 

SSM I Rci = 2.73 ± 0.38SNU, (7) 

SSM II Rci = 2.24 ± 0.23SNU (8) 

we see that both SSMs are fully compatible with the data. That is not the case when 
using for 612, Am^i the central values of standard LMA, that correspond to a negligible 
magnetic moment, from global solar analysis and KamLAND [HI [3T] : 

9 = 0.592, Am^ = 8.0 x 10"^e1/l (9) 

In fact, in this case the Chlorine rate is fully compatible with the SSM II model while 
SSM I is slightly disfavoured at 1.3a. The predictions obtained from the solar neutrino 
and KamLAND global analysis central values are shown in Table ([3]). 





R (SNU) 




SSM I 


SSM II 


pep 


0.120 ±0.002 


0.123 ±0.001 


'Be 


0.638 ±0.067 


0.572 ±0.053 


'B 


2.351 ±0.376 


1.864 ±0.224 


13iV 


0.028 ±^,Z 
0.082 ±0:OiI 


Q Qiq -uo.003 

U.UiJ =1=0.008 
n QK] -UO.002 
U.UOi =1=0.007 


J2 Rci 


3.22 ±0.45 


2.63 ±0.27 



Table 3: Expected event rates in Chlorine for standard LMA with Solar±KamLAND 
bestfit parameters: 9 = 0.592, Am^ = 8.0 x 10"^. 
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3 The SNO+ experiment 



Among the existing and proposed solar neutrino experiments to come online in the near 
future, SNO+ will be the only one with the ability to measure a survival probability, as 
the ratio between the measured and SSM predicted rate of a solar neutrino component, 
at the precision level of 5%. Due to the depth of SNOLAB, SNO+ will have a low 
level of ^^C background, allowing for an accurate measurement of the the pep neutrino 
flux - predicted in the SSM with an error of (1-1.5)%. This is an advantage over the 
measurement of the ^Be flux, that can be done earlier in other experiments (Borexino, 
KamLAND), but is predicted with an error of (9.3-10.5)%. In addition, SNO+ will also be 
the best experiment to measure the CNO-cycle fluxes. Measurements of the pep and CNO 
fluxes will be possible also in Borexino [38]. although with larger systematic and statistical 
uncertainties than at SNO+. The expected accuracy on these fluxes will allow for the 
distinction between different Solar Models accepted by the present data. In addition, 
SNO+ will measure reactor anti-neutrinos, geo-neutrinos and can later be upgraded to 
detect neutrinoless double-beta decays, in order to search for the absolute neutrino mass. 

3.1 The SNO+ detector 

SNO+ [32l |33] is a proposed upgrade of the the Sudbury Neutrino Observatory (SNO) 
detector [31], in which the heavy water target will be replaced by an organic liquid scintil- 
lator. The scintillator was chosen to have a good light yield and transparency, and to be 
compatible with the existing SNO components: chemical and optical compatibility with 
the acrylic vessel and emission wavelength peaks close to the PMT response. 

SNO-I- has a fiducial volume of 1000 tonnes, in a 12m diameter acrylic vessel viewed by 
9456 PMTs mounted on an 18 m diameter geodesic structure. The PMTs have a diameter 
of 20 cm, and are coupled to optical reflectors, increasing the total effective coverage to 
54%. The region outside the acrylic vessel is filled with 7000 tonnes of light water (1700 
tonnes within the PMT structure), acting as a shield against external radiation. The 
external regions viewed by additional 114 PMTs, to act as a veto for cosmic ray muons. 

The large volume of the detector will allow for an effective veto of external backgrounds 
(mainly gammas and neutrons) through position reconstruction, while maintaining a large 
fiducial volume. Concerning internal radioactivity, the purity levels achievable for scintil- 
lators can be estimated from the existing KamLAND ones for the U and Th chains, while 
others such as "^^K and "^^^Bi are being studied for several experiments and expected to 
be reduced by several orders of magnitude. 

In addition, the center of the SNO detector is at a depth of 2092m, or 6010 m of 
water equivalent - it is located in the deepest underground physics facility, SNOLAB. The 
^^C contamination produced by cosmic ray muons - which in general prevents or severely 
hinders the measurement of the pep solar neutrino line - is not a problem at SNO+ since 
at this depth there are only approximately 70 muons to enter the detector per day, and 
so the data taken just after can be cut away. 

The fact that SNO+ is located in a laboratory whose background conditions are well 



7 



known and has the same geometry as the SNO detector, allows for accurate estimations 
of its capabilities even before construction. 

3.2 Neutrino Measurements at SNO+ 

Neutrinos interact in the scintillator through elastic scattering off electrons (z/e~ ve~). 
In the energy range of pep solar neutrinos, the cross-section is around five times smaller for 
muon or tau neutrinos than for electron neutrinos, so SNO+ is primarily a "disappearance 
experiment" for electron neutrinos, even if there is some sensitivity to the other flavors of 
active neutrinos. 

The high light yield of the scintillator - about 100 times more light than the Cerenkov 
light in heavy water - allows the detection of recoil electrons with energies as low as tens 
of keV. 

Since the measured energy will be roughly proportional to the number of detected 
scintillation photons, the main contribution to the energy resolution can be written as: 




where Nph is the detected light yield which, from conservative Monte Carlo simulations 
[35] . was estimated at about Nph = 600 photons per MeV. This results in an energy 
resolution of 4% at 1 MeV, better than in other large liquid scintillator detectors - 6.2% 
at Kami AND [7j and 5% at Borexino[36]. 

The energy of the incoming neutrino is not directly reconstructed, but the different 
electron spectra structures - and mainly the Compton edge directly related to the neutrino 
energy - can be used to statistically separate the neutrino signal from background and the 
different contributions to the solar neutrino flux, namely the mono-energetic pep line. 

The high rate low energy background from ^^C in the scintillator and the gamma 
background from detector materials can be suppressed with data selection cuts in both 
energy - threshold of 500 keV - and position. For the pep analysis, an energy window 
from 800 keV to 1500 keV is expected, to avoid background from fluctuations of the ^ Be 
signal on the lower side. 

Assuming U and Th background levels from KamLAND and the expected levels for 
'^^K and '^^'^Bi after purification, SSM shapes for CNO neutrinos and LMA oscillation, 
a likelihood fit to the energy spectrum allows the separate measurement of the number 
of pep (CNO) neutrino events with a 4% (6%) uncertainty [35], after three years of data 
taking. Adding in a global systematic uncertainty on the fiducial volume, estimated at 
3%, a total flux measurement error of 5% is obtained for pep. 

4 Sensitivity of SNO+ to RSFP and solar models 

In this section we calculate the predicted event rates of pep and CNO solar neutrinos in 
SNO+ for the LMA+RSFP model in the high and low field cases, for the two sets of solar 
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models, SSM I and SSM II (Sec. I4.ip . We then use the predicted rates to estimate how 
sensitive is SNO+ to an RSFP effect (Section 14.21) . and to the discrimination between the 
two solar models, SSM I and SSM II (Section l4.3p . 




Figure 2: Standard LMA and LMA+RSFP spectra in the active sun period for type I 
model (SSM I [23]). 



4.1 Expected rates and spectral analysis 

The expected event rate for each solar neutrino component can be given by 

R{3) = VfidPe dX^ dTF{T,T^) dE <P,^,{E) P.^E) "j^' ^ (11) 

The total fluxes 0j for each component j G {pep, A^, 0} are taken from table ([1]). ^j{E) 
is the normalized neutrino flux spectrum for flux j (delta function for pep). Eq. (ITT!) 
includes a sum over x G {e, /i, r}. The quantity Tm is the observed electron kinetic energy. 
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Figure 3: Same as fig. 3 for type II model (SSM II [21] )• 

T is the true electron kinetic energy, E is the neutrino energy with Emax given by the 
kinematics and Tmin, Tmax are the lower and upper thresholds of the Tm analysis energy 
window, given in Section [31 pe is the electron density of the scintillator and Vfid is the 
analysis fiducial volume. Py^{E) is the probability for z/g conversion, day^{E,T)/dT 

is the cross section for scattering and -F(T, T^) is the response function of the detector. 
This is given by 

where a is the detector energy resolution, as described by Eq. [TOl For models I and II we 
compute the contribution from each component of the solar neutrino spectrum: pep, ^^O, 
and ^^A^ (the F contribution is negligible and is not evaluated). 
We perform the rate calculation for four scenarios: 

SFP The oscillation probabilities Pi,^ are computed with the RSFP effect in addition to 
LMA as explained in Section 12.11 and in [TSj [16] . We consider a high field value of 
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Figure 4: Standard LMA spectrum for models of type I (SSM I [23]) and type II (SSM II 

m)- 

Bq = 280kG, which would correspond to the first three years of data taking, since 
SNO+ is expected to start during the next period of rising solar activity including 
its maximum around 2011 [|3 7]. 

SFPO Pi,^ are obtained from the LMA oscillation effects with Am^, 6 as in eq. 
corresponding to Bq = OkG, the RSFP prediction for low solar activity periods. 

LMA P^^ are obtained from the LMA oscillation effects with Solar+KamLAND bestfit 
Am^, 9, as in eq. Qj. 

NOsc For reference we also evaluated the case for the absence of oscillations, where 
P = S 

The extraction of the pep and CNO signals from the future experimental data will 
be based on a fit to the measured energy spectrum, which will require a very accurate 
knowledge of the detector response, as well as of the residual backgrounds that pass the 
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selection cuts, based on extensive detector calibrations. The SNO+ collaboration has 
carried Monte Carlo simulations [35] of the expected backgrounds in the pep-CNO energy 
window, and performed maximum likelihood signal extraction on the simulated signals and 
backgrounds, mostly from isotopes of the '^^^JJ and '^^^Th, as well as ^^K. The resulting 
sensitivity strongly depends on the background levels. Assuming the target values for 
KamLAND, a sensitivity of 4% for pep and 6% for CNO is obtained for three years of 
data. 

This calculation assumed standard LMA oscillations, and might be changed in the 
case of RSFP, in which the signals are reduced. We carry out a simple spectral analysis, 
without considering the backgrounds, for all cases. 

The energy spectrum, or expected number of events in the spectral bin i is given by: 

Ni = e dTm dTF{T,Tm) dE <Pj<^,{E) P,^{E) ^ (13) 

JTi Jo Jo dl 

where e = Livetime x VjidPe is the exposure, and the sum extends over the flux index 
j e {pep,N,0}. 

The contributions to uncertainty {5Ni) include: 

• the statistical error, considering the number of events N for 3 years; 

• the energy scale error (we assumed an error of IQkeV in the determination of the 
threshold) ; 

• a global systematic error of 3% from the fiducial volume determination; 

• the error in the total flux predictions as given in table ([1]). 

The error for each bin is obtained by adding quadratically these four sources of error: 

5N, = ^(v^)' + (5iVfcale)' + (jy, * 3%)' + (6N^'^^)\ (14) 

where the first three terms are experimental errors (statistical, energy scale, and global 
systematic uncertainty) and the last is the theoretical one (flux uncertainty). Our ex- 
traction uncertainties are in reasonable agreement with the values quoted in [35], before 
including the theoretical and energy scale errors, and can thus be taken as indicative val- 
ues. As expected, the CNO extraction uncertainty is increased in the RSFP case, and that 
is taken into account by this procedure. Equation [13] can be rewritten in matrix form as: 

A^^ = 0,. (15) 

Here the normalized predicted spectra of each component Fij is a 100 x 3 matrix. Once 
the data on Ni are known, one can extract <pj by inverting F: 

0, = Fjl' (16) 
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where the 3 x 100 matrix F ^ = (^F'^ Fj F'^ is the pseudoinverse of F and the transpose 
of F, F'^, is a 3 X 100 matrix. The errors are also calculated from this matrix inversion, 
assuming full correlation (uncorrelation) between the systematic (statistic) errors in each 
bin. 

The spectra used in this analysis are shown in Figures [2}|31 The rate results for type I 
and II models are shown in Tables ([4]) and fl5|) . 





Number 


of events (in thousands) from SSM I 


Component 


RSFP 


RSFPO 


LMA 


NOsc 


pep 


4.82 ±0.28 


6.55 ±0.38 


6.09 ±0.36 


9.98 ±0.58 




(±0.11) 


(±0.15) 


(±0.14) 


(±0.23) 


13 JV 


0.26 ±0.09 


0.40 ±0.14 


0.38 ±0.13 


0.60 ±0.21 




(±0.03) 


(±0.05) 


(±0.05) 


(±0.07) 


15Q 


2.61 ±0.86 


3.63 ± 1.20 


3.38 ± 1.11 


5.52 ± 1.82 




(±0.09) 


(±0.13) 


(±0.12) 


(±0.20) 


Total 


7.69 ± 1.04 


10.58 ±1.45 


9.85 ±1.35 


16.09 ±2.20 




(±0.23) 


(±0.33) 


(±0.31) 


(±0.50) 



Table 4: Expected number of events in 3 years of SNO± from pep, ^^O and ^^N solar neutrinos 
for three years of data taking, assuming the type I model, considering LMA±RSFP, at high 
field (RSFP) and at low field (RSFPO), standard LMA and the no oscillation case. The error 
in parentheses is from the energy scale uncertainty only. The total error is indicated next 
to the predicted value and includes also the uncertainties from the model, from statistics and 
systematics, as stated in the text. 





Number of events (in thousands) from SSM II 


Component 


RSFP 


RSFPO 


LMA 


NOsc 


pep 


4.92 ±0.29 


6.68 ±0.38 


6.22 ±0.36 


10.19 ±0.58 




(±0.11) 


(±0.16) 


(±0.15) 


(±0.24) 




0.17 ±0.03 


0.26 ±0.05 


0.25 ±0.05 


0.39 ±0.08 




(±0.02) 


(±0.03) 


(±0.03) 


(±0.05) 




1.63 ±0.27 


2.26 ±0.38 


2.11 ±0.35 


3.44 ±0.57 




(±0.06) 


(±0.08) 


(±0.08) 


(±0.12) 


Total 


6.72 ± 0.48 


9.25 ±0.67 


8.58 ±0.62 


14.02 ± 1.01 




(±0.19) 


(±0.27) 


(±0.25) 


(±0.41) 



Table 5: Same as Table H for SSM II. 
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4.2 Analysis of the pep results 

The most relevant aspect of the pep resuhs is the large difference between the LMA and 
RSFP predictions, that is not significantly affected by our present uncertainty on the solar 
model. 

In fact, even if the the lowest central value for the best-fit LMA prediction (in SSM I) 
is measured, it is high enough to exclude the highest RSFP value (in SSM II), at more 
than 4(7, and if it would be down by la of the prediction, RSFP would still be disfavoured 
at around 3a. 

On the other hand, if SNO+ finds a pep flux lower than 1.40 x 10^ cm~^s~^, it will 
necessarily follow that new physics is at work. If the RSFP effect is observed, and the 
central value predicted in SSM II for Bq = 280kG is measured, the best fit LMA point is 
excluded at 3.3a. In this case, all the allowed 90% C.L. region of LMA is also excluded 
at the 3a level. 

In addition to a low pep measurement, RSFP predicts significant time variations of 
the measured flux with a solar cycle periodicity, due to the dependence of the survival 
probability on the magnetic field peak. This expected variation would allow the distinction 
between RSFP and other scenarios that predict lower rates regardless of the magnetic field. 



The effect is shown for field values different from Bq = 280kG in figure [5], for A 



01 



values around the resonant value obtained from [291 SO]. We plot in fig [5] the expected 
rate reduction for the pep flux in relation to the non-oscillation function of the 

solar magnetic field and Attiq^. From tables (j4]) and ([5]) it is readily seen that Bq = (the 
X-axis in fig [5]) corresponds to a rate reduction of 65%. 



3.x 10-^ 2.x 10-^ 



1.x 10-^ 



300K 



I 200K 
O 

m lOOK 




55% 
^60% 



300K 



200K 



lOOK 



3.x 10 



-8 



2.x 10 



,-8 



l.xlO"^ 



Am§i (eV^) 



Figure 5: The expected rate reduction for the pep flux with respect to the non-oscillation 
case, as a function of the peak value Bq of the solar magnetic field and AttIq^. 



14 



4.3 Analysis of the CNO results 

RSFP would change also the CNO rates in a consistent way, and that could serve as an 
independent cross-check of the oscillation model. However, the prediction for CNO rates 
depends strongly on the SSM considered and has a large uncertainty within SSM 1. 

In fact examining, for example, the ^^O contribution in the LMA case and taking 
into account both the theoretical and experimental errors, we see that model I predicts 
3383 ± 1047 events, while model II predicts of 2109 ± 361 events, so that the two results 
overlap only slightly. This is also the case for the RSFP predictions, so that after the 
oscillation pattern is established with pep, the distinction between SSMs is improved. 

If the central values of SSM I would be measured, SSM II could be excluded at around 
2.6(7, both in LMA and RSFP cases. The la upper edge of the prediction would allow 
a full exclusion at more than 4. 60", while the lower edge is compatible with SSM II. The 
large uncertainty in SSM I makes is harder to exclude, even if the measurements are in 
full agreement with SSM II predictions. The results could, in turn, be used to constrain 
the upper values of CNO fluxes allowed in this model. 

We see that SNO-I- will be able to determine the ^^N and ^^O fluxes with such a 
precision that it can discriminate SSM of type I from SSM of type II, or severely constrain 
SSM I. There is no doubt that the discrimination of the several flux components with this 
precision will have an impact on our present knowledge of the solar models. However, an 
answer to the question How well SNO+ discriminates SSM I from SSM 11 can be obtained 
only after the experimental data are available. 

4.4 SNO+, Borexino and KamLAND 

We conclude this section with a brief comparison of the model predictions for these three 
real time liquid scintillator experiments. We used the parameter values considered before 
(Am^o = 1.7 X IQ-^eV^, Bq = 280kG), the field profile as in eqs. (2), (3) and Am^i, 9 as 
in eq.(4). The results for the ratio between the predictions of RSFP and LMA with low 
field are shown in table (6) for Borexino and SNO-I- (central values only). For KamLAND 
the same ratio for Be gives 66.0% (see ref.p6j). 

We observe a sizable magnitude difference between Rse and Rpep- This is due to the 
fact that the Be and pep neutrino lines resonate at 0.68 and 0.78 of the solar radius where 
the field is respectively 99% and 61% of its peak value. All other results are consistent 
between Borexino and SNO-I-, taking the errors into account. 

5 Conclusions 

One of the key questions that the SNO-I- experiment will be able to address is the distinc- 
tion between the two classes of SSMs, type I and II. Moreover SNO-I- has the potential 
to distinguish the important issue of variability of the solar neutrino flux for low and 
intermediate energies. In order to investigate the prospects for the SSM distinction we 
looked for an event rate prediction whose model dependence leads to two well separated 
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Component 


Borexino 


SNO+ 




65.7% 




Rpep 


70.9% 


73.6% 


Rn 


66.7% 


65.0% 


Rq 


67.0% 


71.9% 



Table 6: Ratios between RSFP and LMA event rate predictions for Borexino (evaluated as in 
[13]) and for SNO+ (from tables 4 and 5). 

answers. The total flux measured by SNO stands in between the SSM predictions 
and so cannot resolve the ambiguity. Furthermore the Chlorine data can also not pro- 
vide a clear distinction between the two models, the main reason being that the Chlorine 
data combine several intermediate and high energy fluxes (^Be, CNO and '^B). In fact for 
SSM I, while the LMA+RSFP prediction is fully compatible with the data, the standard 
LMA one is disfavoured at 1.3a. On the other hand for SSM II both scenarios (LMA and 
LMA+RSFP) are equally consistent with experimental evidence. 

SNO+ will be able to accurately measure the pep and CNO fluxes. The former, largely 
independent of solar models, will supply the survival probability at low energies, essential 
to distinguish standard LMA from LMA+RSFP. Consequently SNO+ will be able to 
severely constrain the RSFP interpretation, thus strongly favouring LMA or vice-versa. 
The CNO measurement will on the other hand favour one SSM with respect to the other. 
Thus four cases can be classifled according to whether LMA+RSFP or standard LMA is 
favoured by the pep measurement and SSM I or SSM II is favoured by the CNO one. 
We have seen that, if Amg^ is close to 1.7 x 10^^ eV, the value corresponding to the 
most efficient resonances, SNO+ will not only be able to discriminate standard LMA from 
LMA+RSFP after three years of data taking, but might also discriminate SSM II from 
SSM I. 

Would the RSFP explanation be ruled out - or severely constrained - by the pep mea- 
surement, and the LMA interpretation favoured, the Chlorine results favour the SSM II 
solar model, with low heavy element abundance. The CNO results could then be used to 
further confirm the consistency of the model. If, on the other hand, RSFP is favoured 
by the SNO+ pep-data, CNO can then be used to identify the right solar model - in- 
distinguishable from Chlorine and high energy neutrino fiux data alone. For the RSFP 
case considered, if the SSM I prediction for CNO fiuxes is found, the SSM II model can 
be excluded at 3.5a. Although it will be harder to discriminate SSM I from SSM II, the 
allowed regions of SSM I could be severely constrained. In any case, other measurements 
of low energy solar neutrino rates could complement the identification of the solar model 
and be used to test its consistency, and the consistency of the favoured oscillation pattern. 

In particular, consistency between the future pep and CNO measurements at SNO+ 
and the Chlorine experiment requires that SNO+ will either favour an LMA+RSFP sce- 
nario or indicate a preference for a low CNO fiux prediction (as in SSM II). 
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